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ABSTRACT 

In the standard model of disk galaxy formation, a dark matter disk forms as massive satellites are 
preferentially dragged into the disk plane and dissolve. Here, we show the importance of the dark 
disk for direct dark matter detection. The low velocity of the dark disk with respect to the Earth 
enhances detection rates at low recoil energy. For weakly interacting massive particle (WIMP) masses 
-WwiMP ^ 50GeV/c^, the detection rate increases by up to a factor of 3 in the 5 — 20keV recoil 
energy range. Comparing this with rates at higher energy is sensitive to Mwimp, providing stronger 
mass constraints particularly for AfwiMP ^ 100 GeV/c^. The annual modulation signal is significantly 
boosted and the modulation phase is shifted by ~ 3 weeks relative to the dark halo. The variation 
of the observed phase with recoil energy determines Mwimp, once the dark disk properties are fixed 
by future astronomical surveys. The constraints on the WIMP interaction cross section from current 
experiments improve by factors of 1.4 to 3.5 when a typical contribution from the dark disk is included. 
Subject headings: dark matter ~ cosmology. Galaxy: formation 



6 



> 

o\ 

00 
(N 

^' 

O 

oo 

o 



X 



1. THE MILKY WAY'S DARK DISK 

A mysterious dark matter dominates the matter con- 
tent of the universe. Although there are no dark mat- 
ter candidates in the standard model, they are plenti- 
ful in extended models. Among these, weakly interact- 
ing massive parti c les (WIMPs: Cee fc Weinberg Il977t 
iGunn et al\ 119781: lEllis et al\ [198j)), which may arise 
in supers ymmetric extensions of the standard rn odel 
(SUSY) (jJungman. Kamionkowski fc Griesll Il996f l or 
in theories that include universal extra dimensions 
(jChcng. F eng and MatchevI 120021 : iHooper and Profumol 
[2007), stand out as well motivated and detectable. 

WIMPs may be d etected directly by scatterin g in a 
laboratory detector (| Goodman and WittenI [l985l ) or in- 
directly by their an nihilation products from their high - 
est density regions (jSilk fc Srednickil[l98l lLak3fl99l . 
In the case of direct detection, we must know the dark 
matter's phase space structure to predict rates. In early 
calculations, the standard halo model (SHM) of the dark 
matter assumed no rotation and the density distribu- 
tion was taken to be a spherical isothermal sphere with 
a core radius of several kpc. More recent modeling in- 
clude s the cuspier central profiles from ACDM simula - 
tions avarro. Frenk fc Whit3IT997HMoore et aZ.|[l998l) . 
producing changes of 0(10%) with re spect to the SHM 
([Kamionkowski fc K inkhab walal Il998f ). Larger boosts 
have been c laimed if dark matter is highly clumped 
(lGreenll2002D . but it is more likely that we live outside a 
clump, leading to a modest reduction in th e local density 
([c. f. Kamionkowski fc Koushiappas Il2008f ). 

Simulations containing only dar matter particles have 
extremely high resolution, but they may not be address- 
ing the "next to leading order" of the model because 
they do not include the effect of the baryons. Read 
et al. (2008) recently demonstrated that massive satel- 
lites are preferentially dragged into the baryonic disk 
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plane by dynamical fricti on where they dissolve l eaving 
a thick dark matter disk ([Read et aLll2008t lLakel [r989). 
The precise properties of the dark disk depend on the 
stochastic merger history and cosmology. However, given 
the expected merger history for a typical Milky Way in 
ACDM, they found a dark disk with density in the range 
Pd/ Ph ^0.2 — 1 at the solar neighborhood (where ph is 
the density of the SHM). The lower bound Pd/ph = 0.2 
is particularly conservative since it is produced by just 
one merger of Large Magellanic Cloud mass within 20° 
of the disk plane. In ACDM we expect two such low in- 
clination mergers per Milky Way (and seven in total at 
all inclinations). 

We may obtain an upper bound on Pd/ph fror n the 
kinematics of stars at the solar neighborhood ([OortI 
11932.19601: [Bahcal]|[T98l . The latest measurements from 
Hipparcos give a conservativ e upper bound of Pd/Ph < 3 , 
including systematic errors (iHolmberg and Fli2m| |2003; 
[Statlcr 1989); with pd/ph being < 2 likely. If more than 
half of the thick disk owes to accretion, the likely dark 
disk density would be near the upper limit. As such we 
consider Pd/ph — [0.5,1,2] in this paper. The disk den- 
sity Pd is an excess over ph, locally increasing the dark 
matter density. 

The kinematic properties of the dark disk can be esti- 
mated from the accreted stellar thick disk that forms con- 
currently. Stellar thick disks are found in the Milky Way 
and in all well-observed spiral galaxies fBurstein' '1979|; 
iGilmore and ReidI Il983t [^achim fc Dalcanton 200§,), 
while at least one counterrotati nq thick disk is strong 
evide nce for an accretion origin ([Yoachim fc Dalcanton I 
|2005| ). However, thick disks ca n also form through heat - 
ing of an underlying thin disk ([Kazan tzidis et a/."2007'). 
or even directly from extended gas (Brook et al. 2004); 
indeed it is difficult in ACDM to o btain a thick disk mas- 
sive enough from accretion alone ([Read et aZ.l i 20081. In 
this paper, we assume - based on the numerical mod- 
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els of |Read et al\ I2008D - that the dark disk's kine- 
matics match the Milky Way's stellar thick disk. At 
the solar neighborhood, this gives a rotation lag viag 
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Fig. 1. — DifTerontial recoil rates for Gc (red/gray) and Xe 
(blue/black) targets, for AfwiMP = 100 GeV/c^ and o-(--yyii^p = 
10~* pb in the SHM (solid line) and the dark disk. Three different 
values of Pd/Ph (0-5 dashed, 1 X and 2 A) are shown. The vertical 
lines mark current experiment thresholds: XENONIO (blue/black) 
using a Xe and CDMS-II (red/gray) using a Ge target. 
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Fig. 2. — Recoil energy below which the signal is dominated by 
the dark disk (compared with the SHM) as a function of AfwiMP 
for Ge (red/gray) and Xe (blue/black) targets. Three different val- 
ues of Pd/Ph (0.5 dashed, 1 X and 2 A) are shown. The horizontal 
lines mark current experiment thresholds: XENONIO (blue/black) 
using a Xe and CDMS-II (red/gray) using a Ge target. 



of 40 — 50 km s with respect to the local circular ve- 
locity, and dispersions of (iTfl, cr^, CTz) — (63, 39, 39) km 
s-i (|Read et aLll2008D . Since the dispersions are nearly 
isotropic and somewhat uncertain, we model the dark 
disk with a simple one-dimensional Maxwellian distribu- 
tion, with a dispersion and lag, a = viag = 50 km s~^ 
to show its general effect on direct detection. Improving 
on this assumption will involve untangling heated ver- 
sus accreted components in the Milky Way stellar thick 
disk. This should become possible with future astronom- 
ical surveys like the Ra dial Velocity Experiment (RAVE) 
(jSteinmetz et al\l200(^ a nd Global Astrometric Interfer- 
ometer for Astrophysics (jPerrvman e t a/."2001') that will 
provide full 6D phase space and chemical information for 
hundreds of thousands of individual stars. 

We emphasize that the dark disk must form in any hi- 
erarchical model of structu re formation; it is not special 
to our ACDM simulations ("Rea d et alll2008[ ). However, 
ACDM is specific enough to predict the dark disk density 
with uncertainties owing only to the stochastic nature of 
the merging and accretion history. T he dark disk is very 
different from dark matter strea ms (|Freese et al\ l2004t 
ISavage. Freese fc Gondololl2006f l that have a low filling 
factor (we are not likely to live in a stream), and are 
stochastic micro structure. By contrast, the dark disk is 
the expected equilibrium end state of dissolving satellites 
and the Earth must be embedded in one (if hierarchical 
formation is correct). Like the near-spherical dark mat- 
ter halo, the dark disk is macrostructure. 

2. DIRECT DETECTION AND THE DARK DISK 

Direct detection experiments measure nuclea r recoil 
rates above the detector's energy threshold (jBau_disl 
[2006); here we consider Ge and Xe. The detected elastic 
WIMP-nucleon recoils will range from a few to tens of 
keV. The expected recoil rate p er unit mass, unit nuc lear 
recoil energy and unit time is (jLewin fc Smithlfl99l 
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where p is the local dark matter density {ph = 0.3 
GeV cm-3 in the SHM), cr(wiMP,N) is the WIMP- 
nucleus scattering cross section, F{E) is the nuclear 
form factor, Mwimp and m are the masses of the dark 
matter particle and of the target nucleus, respectively, 
fj, is the reduced mass of the WIMP-nucleus system, 
V — |v| and Vmax is the maximal velocity in the earth 
frame for particles movin g at the ga l actic escape ve- 
locity Vesc = 544 km s"^ ijSmith et aLll2007D . We only 
consider the spin-independent (SI) scalar WIMP-nucleus 
coupling in this paper, since it dominates the inter- 
action (depending however on the dark matter parti- 
cle) for target media with nucleon num ber ^ > 30 
pungman. Kamionkowski fc GriesH |1996[ ). We model 
the velocity distributions of particles in the dark disk and 
the SHM with a simple one-dimensional Maxwellian: 



/(v, t) cx exp 
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where v is the laboratory velocity of the dark matter 
particle and the instantaneous streaming velocity — 
Vrirn + V0 -|- Vorb(i)- Thls Streaming velocity is the sum 



of local circular velocity Vcirc 
p ecu liar motion of the Sun 



(0,220,0) km 



the 



(10.0,5.25,7.17) km 



(iDehnen and Binnevlll999l ) with respect to Vdrc and 
the orbital velocity of the Earth around the Sun Vorb {t) ■ 
In the SHM, the halo has no rotation and the dispersion 
o" = |vciic|/V2- For the dark disk, the velocity lag viag ~ 
(0, 50, 0) km s~^ replaces Vcirc and a dispersion of 50 km 
s~^ is adopted. 

The lower relative velocities of the dark disk signifi- 
cantly increase the differential rate at low energies com- 
pared with the SHM rate (Figure 1). Detection of the 
dark disk crucially depends on the detector's low energy 
threshold. The differential rate for a specific WIMP tar- 
get depends on Mwimp- In Figure 2, we show the energy 
below which the dark disk dominates the rate as a func- 
tion of AfwiMPj for three values of pd/ Ph- The total rate 
in a detector is the sum of the two contributions from 
the SHM and the dark disk, which dominate at high and 
low energies, respectively. For Mwimp ^ 50GeV/c2, the 
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Fig. 3. — Annual modulation shown as the residual counting rate vs. date for the XENONIO experiment (4.5 to 27 keV). The residuals 
are calculated with respect to the mean counting rates (given as numbers over each line) using a WIMP-nucleon cross section of 10~* 
pb. The top/bottom row is calculated for Pd/Ph = 0.5/1 and Mwimp (left to right) of 50GeV/c2, 200GeV/c2 and 500GeV/c2. The 
(blue/dashed) lino is the modulation signal obtained from the SHM, the (red/dot-dashed) line is the modulation signal from the dark disk 
and the (black/solid) line is the total modulation signal. The maximum of the dark disk contribution is shifted to May 9th compared with 
the SHMs maximum/minimum on May 30th. Note the different vertical scales in each of the three columns of the plot array. 
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Fig. 4. — Total differential rate for a Ge target in the SHM (solid) 
and four different values of pd/ph (0-1 dashed, 0.2 X, 0.5 A, 1 
V and 2 □) are shown for two Mwimp (lOOGeV/c^ (red/gray) 
and 200 GeV/c^ (blue/black)). The vertical lines mark the current 
CDMS-II threshold and a threshold of 1 keV. 
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Fig. 5. — Phase shifts as a function of A/wiMP in the energy 
range reported by the CDMS-II experiment (10-100 keV; red/grey) 
and the XENONIO experiment (4.5-27 keV; blue/black), for three 
different values of Pd/Ph (0.5 dashed, 1 X and 2 A). 



dark disk contribution lies above current detector thresh- 
olds, giving a much greater change in detection rate with 
recoil energy compared with the SHM alone. 

The total rate in a detector using a Ge target is shown 
in Figure 4 varying pd / Ph and Af-yviMP • If the detectors 
threshold is sufficiently low even an extremely conserva- 
tive dark disk with Pd/ph = 0.1 can be detected. Cur- 
rent germaniuni detectors achieve thr esholds below 1 keV 
()Lin aZ.ll2008HAalseth et aZ.ll2008[) . The details of the 
differential rate with energy, as shown in Figure 4, be- 
tray both the contribution of the dark disk relative to the 
SHM and AfwiMP- This introduces a mass-dependend 
characteristic shape of the differential rate which will im- 



prove the constrains on A/wimp upon detection. 

The motion of the Earth around the Sun gives rise to 
an annual modulation of the event rate and re coil energy 
spectrum (jDrukier. Freese and SpergeTlllQSGf ). The an- 
nual modulation is more pronounced for the dark disk, 
since the relative change to the mean streaming velocity 
owing to the Earth's motion is larger (^19%) compared 
with the SHM (~6%). We show in Figure 3 the residual 
integrated rates for a liquid xenon detector throughout a 
year, for three different AfwiMP and two values of pd/ Ph- 
The residuals are calculated with respect to the mean 
counting rates in a given energy region. 

The phase (defined at maximum rate) of the dark 
disk and the SHM differ because the Sun's motion is 
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I Angle et aLll2008[ ). Depending on pd/ph, we exclude new 
regions in the allowed parameter space for A/wimp ^ 
50GeV/c2. 

On a final note, we find that including the dark 
disk component does not change the interpretation of 
the annual modulatio n signal observed in the DAMA 
(|Bernabei et aLl [20081) experiment for pure SI coupling. 
At high MwiMP, the allowed DAM A region is still ex- 
cluded by C DMS-II and XENONIO (I Ahmed et aLI 120081: 
lAngle etal . 2008) results, while at low Af-yviMP no new 
parameter region opens. 
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Fig. 6. — Effect of the increased dark matter flux on the SI 
WIMP-nucleon cross section constraints obtained by the CDMS- 
II (red/gray) and XENONIO (blue/black) experiments, for three 
different values of pd/Ph (0-5 dashed, 1 X and 2 A). The solid 
lines give the constraints if only the SHM component is considered. 



slightly misaligned to the dark disk. While the phase 
of each component does not depend on Mwimp, their 
sum does because their amplitudes depend on Mwimp- 
We show this dependency in Figure 5, for three values 
of Pd/Ph- The phase shift is determined by the rela- 
tive contributions of each component. Figure 3 shows 
that the phase shift is largest for low pd/ph, since in 
this case the sum preferentially follows the halo modu- 
lation phase, while for higher Pd/ph the disk component 
dominates the modulation phase. This is a new effect 
introduced by the presence of the dark disk that allows 
AfwiMP to be uniquely determined from the phase of the 
modulation signal, for given pd/ph- Note that there is an 
amplitude flip for the SHM that occurs as A/w^jmp is in- 
creased, which is not seen for the dark disk. As Af-yviMP 
is lowered, the "crossing energy" at which the differen- 
tial rates for minimal and maximal WIMP velocity are 
equal shifts to lower energies. For the dark disk, it re- 
mains close to, or below, current thresholds and so the 
amplitude flip is not seen. 

The effect of the dark disk on current upper limits 
on the SI WIMP-nucleon cross secti on is shown in Fig- 
ure 6 for CDMS-II and XENONIO (| Ahmed et aLl[2Q0l 



3. CONCLUSIONS 

In ACDM, a dark matter disk forms from the accretion 
of satellites. We have shown how its low velocity with 
respect to the Earth alters the expected rate and annual 
modulation signal in dark matter detectors. Our main 
findings are as follows. 

The dark disk boosts the detection rates at low recoil 
energy. For A/wimp ^ 50GeV/c^, recoil energies of 5 
- 20 keV and Pd/ph < 1, the rate is boosted by factors 
up to 2.4 for Ge and 3 for Xe targets. Comparing this 
with the rates at higher energy will constrain A^wimp, 
particularly for A^wimp > lOOGeV/c^. 

The dark disk has a different annual modulation phase 
than the dark halo, while the relative amplitude of the 
two components varies with recoil energy and AfwiMP- 
As a result, there is a new richness in the annual mod- 
ulation signal that varies uniquely with A^Twimp, for 
given dark disk properties (the properties of the dark 
disk will be measured from next generati on surveys 
(jSteinmetz et aLll2006t iPerrvman et aZ.ll2001[ )V 

The increased expected dark matter flux provides new 
constraints on the WIMP cross section from current ex- 
periments. For likely dark disk properties (pd/ph < 1), 
the constraints for pure SI coupling improve by up to a 
factor of 1.4 fo r CDMS-II (lAhmed aLl [200l and 3.5 
for XENONIO (|Angle et aLll2008n . 
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